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Relay communication is considered as a popular solution for expanding 
the coverage, increasing the transmission capacity and reducing the power 
consumption of the communication networks. In this paper, we proposed 
and investigated the two-way decode-and-forward (DF) half-duplex (HD) 
relaying network with the direct link between two sources Si and S2. 
Firstly, the system model, energy harvesting (EH) and information 
transmission (IT) are presented. The closed-form analytical expression 
of the system outage probability (OP) is analyzed and derived in the next stage. 
Finally, the correctness of the analytical expressions is verified by Monte Carlo 
simulation. The research results show that the analytical and simulation 
are the same in connection with the primary system parameters. 
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1. INTRODUCTION 

Nowadays, the fifth-generation (5G) wireless networks is considered as a hot research area 
in academic and industrial based on the increasing data rate transmission. Relay communication, in which 
the relay forwards the signal received by a source to a destination, has a massive consideration in research, 
due to its ability to expand the coverage, increase the capacity, and reduce the power consumption [1-4]. 
Moreover, relaying the communication network is an effective way to combat the performance degradation 
caused by fading, shadowing, path loss, and is an efficient way to improve spectrum efficiency and extend 
coverage [1-6]. On another way, two-way relaying, where two users exchange information with each other via 

a single or multiple relays, provides improved spectral efficiency compared to conventional one-way relaying 
by using either superposition coding or physical layer network coding at relays. Furthermore, to satisfy the 5G 
requirements, relaying schemes with high spectrum efficiency, such as two-way, and full-duplex, etc., have 
been recently attracted considerable attention [7-10]. In [7], the authors investigate one-way full-duplex (FD) 
relaying and two-way half-duplex (HD) relaying to minimize the spectral efficiency loss associated with 
one-way HD relaying. In [8] and [9], the authors proposed and investigated one-way FD relaying with multiple 
antennas to provide a solution to overcome the spectral efficiency loss. In [10], the authors investigate 
one-way FD relaying with opportunistic relay selection to enhance the performance. 
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In this research, the two-way decode-and-forward (DF) HD relaying network with the direct link 
between two sources Si and S 2 is proposed and investigated. The system model, energy harvesting (EH) 
and information transmission (IT) are presented in the first stage. After that, we derived the closed-form 
of the system outage probability (OP) in the next stage. Finally, the correctness of the analytical expressions 
is verified by Monte Carlo simulation in connection with the primary system parameters. The research results 
show that the analytical and simulation are the same in connection with the primary system parameters. 
In this research, we focus on some problems: 

a. The two-way DF HD relaying network with the direct link between two sources S 1 and S 2 is proposed 
and investigated. 

b. The closed-form expression of the system OP is analyzed and derived. 

c. All the results are convinced by Monte Carlo simulation in connection with all primary 
system parameters. 

The structure of the rest of the paper can be drawn as follows. We provide the system model, 
the energy and information transfer phases in section 2. The closed-form expression of system OP analysis 
is derived in section 3. We introduce the results and have some discussions in section 4. In the last section, 
some conclusions are proposed. 


2. RESEARCH METHOD 

In this section, the two-way DF HD relaying network with the direct link between two sources Si and 
S 2 is drawn in Figure 1. The energy harvesting (EH) and information transformation (IT) for this proposed 
model system are proposed in Figure 2. In this protocol, the transmission interval time is T, 
which consists of three-time slots T/3. In the first time slot T/3, the R harvests energy pPi from the source node 

51, and the source uses the energy (l-p)Pi for information transmission to R and S 2 (here 0 < p <\ 
is the power splitting factor). In the second interval time T/3, the R harvests energy pP 2 from the source node 

5 2 , and the source S 2 uses the energy (l-p)P 2 for information transmission to R and Si. Finally, the remaining 
time slot T/3 is used for information transferring from the R to the nodes Si and S 2 [11-16]. 



EH Phase 


IT Phase 



Figure 1. System model 


Figure 2. The energy harvesting (EH) 
and information transmission (IT) processes 


2.1. Energy harvesting phase 

Let Si transmits the symbol xi in the first interval time. We define the received signal at the R node 
R and the S node S 2 as: 

=/ V«T +«' 

Ti,2 — ^1,2T 

Where e||jc,|“| = FJ,E{«} is expectation operator and Pi represents the average transmit power 

at the Si. Further, n' r and n[ denote the zero-mean additive white Gaussian noise (AWGN) with variance No. 
Then, the harvested energy at the R node can be formulated as the following: 
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F 1 



( 2 ) 


In the second interval time, the source node S 2 transmits the signal X 2 to the nodes R and Si. 
Therefore, the received signals at the R and Si can be expressed, respectively, as: 


•E.K }h,R X 2 + n r ’ 

y 2 1 = h 2 X x 2 + n' 1 


(3) 


Where e||jc 2 | 2 | = P 2 , P 2 is the average transmit power at the S 2 . And we assume that n" and n' 1 

are the zero-mean additive white Gaussian noise (AWGN) with variance No. Here, the total harvested energy 
at the R node can be defined as: 


VP 


E h =- 




+ P 2 \h 2J 


(4) 


In this model, the average transmits power from source Si and S 2 is assumed to equal. So, the (4) can be 
rewritten as the following: 


E h 


3 



(5) 


Where we denote P 1 = P 2 = P. Therefore, we have the average transmit power at the R as: 


Pr = 


(773) 


-vpP^Ki 


+ \K 


( 6 ) 


2.2. Information transmission phase 

In the first interval time, after doing EH. Si broadcasts the information to the R node and S 2 with 
remaining power (1 - p)P . From that, we have the received signal at the R node and S 2 as: 


ylR=fi-P\,RXi+K’ 
y’i.2 = Ki* +n 2 


Similar to the first case, the received signal at node and S 1 node can be given as the following: 

3 = ~PE 2 R x 2 +n r , 

(o) 

y‘h = K,i x i + n \ 

Finally, in the third interval time, the received signal at the sources Si and S 2 can be expressed, 
respectively, as: 

•Vl = h R,l X R+ n x ’ 

3^2 — ^R.2 X R ^ 

Where E11w A . | _ J- = P K . The received signal at the source node S 1 or S 2 can be formulated as: 
y? F = ^Vi(*i©*2) + ™ y? F = /PrV;(*i © x r> + n y? F = VV'u.iOi © x 2 ) + n 1 


(9) 


GO) 
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Where i e (1,2). Once again, after applying the canceling the self-interference method and combine with (6) 
the SNR at Si can be defined as: 


n F w (k« r+k*r)ki| 2 =np<m Wi +n ] 


(ID 


Where Vi =\hj 


3. SYSTEM PERFORMANCE ANALYSIS 

Hence, the outage probability (OP) of DF mode at the Si can be formulated as: 

OP DF = Pr {y™ < } = Pr {max {wf>, [<P, + <P t ] ■ ¥<P, } < Y , h } 

= Pr{ww, \v, +<p,]< r lh } xP < r,„} 

Where we denote: 


r 'l 


( 

\ 




Y , 

( 1 n th\ 

| = 1 — exp 


th 

\ W 3 <r rti =ft p < “ j 

V 




3 / 


Which A 3 is mean of the random variable (RV) <p 3 . Here we have: 


( 12 ) 


(13) 


P 3 = Pr{>7 PV<m<Yth\ = Pr|pi< 

f ?th 


Yth l 
nPVV) 


= oH 


np¥<p 


¥\f<p(<p)d<p 


(14) 


Where (p = (p 2 +(p 4 For calculating the probability in (13), the probability density function (PDF) 
and cumulative density function (CDF) of <p is calculated as in Lemma 1 
Lemma 1. 

The CDF of (cp 2 + (p 4 ) can be formulated as: 


F v (x) = Vv[(p<x\ = VY[(p 2 <x-cp 4 ] 


X X-<p 4 X 

= \f v S<P^ d( P4 j f< h ( ( Pi) d( P'l=\ F 9l { X - ( P*)f<pS ( PJ d( P4 

0 0 0 



(15) 


Where A 4 is the mean of RV (p 4 . 


Case 1: A 2 = A 4 = A 
From (28), we have: 


F lp (x) = 1 - exp 



( 16 ) 


Bulletin of Electr Eng & Inf, Vol. 9, No. 1, February 2020 : 265 - 272 





Bulletin of Electr Eng & Inf 


ISSN: 2302-9285 


□ 269 


And the PDF of (p can be computed as: 


f<p(x) 


dF v (x) 

dx 



(17) 


Case 2: A 2 ^ /t 4 

Similar, we can obtain CDF and PDF of <p as followings respectively: 



i 2 

exp 

r 

-exp 

Ml 

l Mj 

M -^2 

l Mj 

l 42 J. 


(18) 


/$>(*)=■ 




e x p |--J- ex p^-- 


We can calculate the probability P 3 into 2 case 
Case 1: A 2 = A 4 = A 

Substituting (14) into (16), we can obtain: 


A = 1 - J exp 


Yth 


= 1 


0 

- —« ) (axexp 

A' 0 


vpwvh 


^exp 


-9 


Tth 


dtp 


|xexp| ~jY<P 


npvvh _ 

Apply (3.471,9) in [17], (20) can be rewritten as: 

As=l— 2 ^—K 2 [2 m 

tipy/A.\ I Vr/pyMi 


(19) 


( 20 ) 


( 21 ) 


Where K v (•) is the modified Bessel function of the second kind and v ,h order. 

Also, finally, substituting (13), (21) into (12), the system OP at the source Si in the first interval time can be 
obtained as: 


OP DF = \ X ~^ ?! 


Yth 


lYth 


hv)\ npv^h 


K 2 \2. 


m 

ripysMl . 


( 22 ) 


Case 2: A 2 ^ A 4 

Substituting (19) into (14), the probability P 3 can be expressed as: 


P, = 1 - J expj 
0 


=1 


Yth 

1 

f 

_<p I 

_ exD fM 

*lPV<Ph_ 

M~^2 

exp 

MJ 

p l 4 2 JJ 


dtp 


—T ! exp 
M-*i 0 


m 


1 


0 

1 


exp 


vpvvh 

Yth 


exp —— -exp —— 


3 3 -f exp 

/I4-/I2 o 


Wf'll 
Tth 

riPWh _ 


/I 4 


xexp —— \dtp 


h 


dtp 


xexp 


a 2 J 


dtp 


Apply eq (3.324,1) in [17], (23) can be rewritten as: 


(23) 
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P}= 1 - 


Yth 


M-^'inpvh 


JMxkJiJ 


m 


Yth 


riPYhh 


(24) 


And then substituting (13), (24) into (12), the OP at the source Si in the second case can 
be obtained as: 



T 

to 

X 

s- 

X 

4Mxk\ 2 / rth ) 

v \ inPYhM J 


[ l %iv)\ 

M~^2 V r lP¥\ 

[TT''Ki f n 1 r,h \ 




A j 



(25) 


4. NUMERICAL RESULTS AND DISCUSSION 

In this section, we investigate the system OP using Monte Carlo simulation in connection with 
the primary system parameters [18-25], In Figure 3, the effect of the energy coefficient r| on the system OP 
is plotted with the primary system parameters as p=0.8, R=0.5 bps/Uz, t|/=0.5 dB. From the results, we can see 
that the OP of Si significantly decreases with the rising of rj. Here, we consider both cases with 
and without a direct link between the sources Si and S 2 . From Figure 3, we can see that the system OP 
with the direct link is better than the system OP without the direct link between the sources. Furthermore, 
the simulation results agree well with the analytical results. In the same way, the system OP versus vp 
is drawn in Figure 4. In this simulation, we set r|=0.5, R=0.5bps/Hz, p=0.5, and consider both cases with 
and without a direct link between the sources. As shown in Figure 4, we can state that the system OP has 
a considerable decrease when vp varies from-10 to 10. In Figure 4, the simulation and analytical results 
are the same with all values ofa and bi, and again, the OP with the direct link is better than the cases without. 
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DF(case 2-presence of direct link) with A 4 =2 and A 2 =1 

s\ 

♦ DFfcase 2-no direct link) 
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* 


-2 0 2 
V (dB) 


Figure 3. OP versus r| 


Figure 4. OP versus vp 


Furthermore, we investigate the effect of R and p on the OP of the model system, as plotted in Figures 
5 and 6. In these Figures, we set the primary parameters, as shown in these Figures. From Figures 5 and 6, 
it can be observed that OP has a slight decrease when p varies from 0 to 1, and the system OP crucially increases 
with the rising of R from 0.5 to 5. In all two Figures, the analytical and simulation results agree well with each 
other. On another hand, the comparison of the system OP of both cases with and without a direct link between 
these sources is demonstrated in all Figures. From the results, we can state that the system OP with the direct 
link is better than another case. 
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OP of system versus R with ?>=0.8,^=0.5 and «/>=5dB 
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Figure 5. OP versus R 


Figure 6. OP versus p 


5. CONCLUSION 

In this research, the two-way DF HD relaying network with the direct link between two sources Si 
and S 2 is proposed and investigated. The system model, EH and IT are presented in the first stage. 
After that, we derived the closed-form expression of the system OP in the next stage. Finally, the correctness 
of the analytical expressions is verified by Monte Carlo simulation in connection with the primary system 
parameters. The research results show that the analytical and simulation are the same as the primary system 
parameters. The results can provide a critical recommendation for the relaying communication network. 
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